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ABSTRACT 

We use the Fifth Data Release of the Sloan Digital Sky Survey to study X-ray-selected galaxy groups 
and compare their properties to clusters. We search for infall patterns around the groups and use these 
to measure group mass profiles to large radii. In previous work, we analyzed infall patterns for an 
X-ray-selected sample of 72 clusters from the ROSAT All-Sky Survey. Here, we extend this approach 
to a sample of systems with smaller X-ray fluxes selected from the 400 Square Degree serendipitous 
survey of clusters and groups in ROSAT pointed observations. We identify 16 groups with SDSS DR5 
spectroscopy, search for infall patterns, and compute mass profiles out to 2-6 /i Mpc from the group 
centers with the caustic technique. No other mass estimation methods are currently available at such 
large radii for these low-mass groups, because the virial estimate requires dynamical equilibrium and 
the gravitational lensing signal is too weak. Despite the small masses of these groups, most display 
recognizable infall patterns. We use caustic and virial mass estimates to measure the scaling relations 
between different obscrvablcs, extending these relations to smaller fluxes and luminosities than many 
previous surveys. Close inspection reveals that three of the groups are sub-clusters in the outskirts of 
larger clusters. A fourth group is apparently undergoing a group-group merger. These four merging 
groups represent the most extreme outliers in the scaling relations. Excluding these groups, we find 
Lx oc <7p 4±1 ' 6 , consistent with previous determinations for both clusters and groups. Understanding 
cluster and group scaling relations is crucial for measuring cosmological parameters from clusters. 
The complex environments of our group sample reinforce the idea that great care must be taken in 
determining the properties of low-mass clusters and groups. 

Subject headings: galaxies: clusters: individual — galaxies: kinematics and dynamics — cosmology: 
observations 



1. INTRODUCTION 

A large fraction of all galaxies are members of groups. 
Groups are a less extreme and much more common type 
of system than galaxy clusters. Because they are less 
massive than clusters, groups are currently less well un- 
derstood than clusters. Here, we investigate the optical 
properties of a sample of X-ray-selected groups to deter- 
mine whether they can be modeled as scaled-down ver- 
sions of clusters 5 . In particular, we study the outskirts of 
the groups to determine if their infall regions are read- 
ily identifiable as they are in the outskirts of clusters. 
With large samples of spectroscopic members, we then 
estimate the virial masses of the groups and determine if 
the groups obey the same scaling relations as clusters. 

Galaxy clusters and groups are surrounded by infall 
regions in which galaxies are bound to the system but 
are not in equilibrium. The Cluster Infall Regions in 
the Sloan Digital Sky Survey (Rines & Diaferio 2006, 
hereafter CIRS) project showed that X-ray-selected clus- 
ters display a characteristic trumpet-shaped pattern in 
radius-redshift phase space diagrams. These patterns, 
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termed caustics, were first predicted for simple spherical 
infall onto clusters (Kaiser 1987; Regos & Geller 1989), 
but later work showed that these patterns reflect the dy- 
namics of the infall region (Diaferio & Geller 1997, here- 
after DG) and (Diaferio 1999, hereafter D99). CIRS and 
earlier similar studies (e.g., Geller et al. 1999; Rines et al. 
2003) showed that the amplitude of the caustics yields 
an estimate of cluster mass profiles consistent with both 
virial and X-ray mass estimates where the techniques 
overlap. More recently, Diaferio et al. (2005) and Lemze 
et al. (2009) showed that caustic masses also agree well 
with mass estimates based on gravitational lensing. Be- 
cause neither galaxies nor gas is expected to be in equi- 
librium outside the virial radius of a cluster, the caus- 
tic technique and weak lensing are the only well-studied 
methods for determining cluster mass profiles at large 
radii (sec Diaferio 2009, for a recent review of the caus- 
tic technique). 

CIRS showed that infall patterns are ubiquitous in 
nearby massive clusters selected by X-ray emission. The 
CIRS clusters are fairly massive clusters and generally 
have little surrounding large-scale structure (but sec 
Rines et al. 2001, 2002). One might suspect that the pres- 
ence of infall patterns is limited to massive, isolated clus- 
ters. However, other investigators have found infall pat- 
terns around the Fornax Cluster (Drinkwater et al. 2001), 
the Shapley Supercluster (Reisenegger et al. 2000), an 
ensemble cluster comprised of poor clusters in the Two 
Degree Field Galaxy Redshift Survey (Biviano & Girardi 
2003), and even the galaxy group associated with NGC 
5846 (Mahdavi et al. 2005). Here, we extend the study 
of infall patterns by studying a large number of systems 
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with smaller X-ray fluxes than the CIRS clusters. Be- 
cause these systems typically have smaller masses and 
fewer member galaxies (more typical of groups), infall 
patterns might not be identifiable in these systems. 

In particular, we use the new 400 Square Degree (400d) 
survey of clusters and groups in ROSAT PSPC pointed 
observations (Burenin et al. 2007). The 400d survey is 
the largest area cluster survey extending to flux limits 
of wl.4x 10~ 13 erg s _1 cm -2 . This flux limit is a factor 
of 20 smaller than the flux limit of catalogs based on 
the ROSAT All-Sky Survey (RASS; Voges et al. 1999). 
We search for 400d clusters and groups in the spectro- 
scopic footprint of the Sloan Digital Sky Survey Data 
Release 5 (Adelman-McCarthy et al. 2007, SDSS DR5). 
This approach is similar to CIRS and to the RASS-SDSS 
(Popesso et al. 2004, 2005) analysis of clusters in DR2. 
Compared to both the CIRS and RASS-SDSS samples, 
the 400d groups are expected to have significantly smaller 
masses on average. 

One motivation for this study is to improve our un- 
derstanding of cluster scaling relations, in particular to 
extend these relations to the regime of high-mass groups. 
Ambitious cluster surveys like the South Pole Telescope 
and Atacama Cosmology Telescope require a good under- 
standing of cluster scaling relations to accurately mea- 
sure cosmological parameters from the abundance and 
evolution of clusters. Quantifying the scatter in these re- 
lations and any Malmquist-like bias is crucial for obtain- 
ing robust cosmological constraints from cluster surveys 
(e.g., Stanek et al. 2006). Many forecasts for dark energy 
constraints from cluster surveys assume that scatter in 
the mass-observable scaling relations has a log-normal 
distribution (e.g., Mantz et al. 2008). Deviations from 
this assumption could significantly impact the constrain- 
ing power of these surveys. 

Recent studies of scaling relations that extend into the 
group regime have often reached contradictory conclu- 
sions. For instance, a common scaling relation is the 
Lx — 0p relation between the X-ray luminosity Lx of a 
cluster or group and its projected velocity dispersion a p . 
For massive clusters, this relation is typically found to 
have a steep slope of 4. 4^3 (Mahdavi & Geller 2001) or 
3.7±0.3 (Popesso et al. 2005). For groups, various stud- 
ies have found slopes as shallow as 0.37 (Mahdavi et al. 
2000, the fit is a broken power-law with a slope of 4.02 
for clusters) and as steep as 4.7±0.9 (Helsdon & Pon- 
man 2000). A detailed optical-X-ray study by Osmond 
& Ponman (2004) recently found a slope of 2.5±0.4. The 
variety in slopes is produced by many factors, including 
differing definitions of Lx, and in some cases, possible 
evolution of the galaxies via dynamical friction (Helsdon 
et al. 2005). One problem is that few large, complete 
group catalogs are available, so many existing studies uti- 
lize either heterogeneous samples or include many groups 
with limited data on each group. One solution is to uti- 
lize X-ray selection (e.g., Mulchaey & Zabludoff 1998; 
Mahdavi et al. 2000), but existing cluster/group catalogs 
based on RASS (Popesso et al. 2005, CIRS) have been 
restricted to systems with relatively high X-ray fluxes 
(and therefore including few groups). The fainter X-ray 
flux limits of the 400d survey allow detection of several 
groups. Combining this catalog with optical data from 
the Sloan Digital Sky Survey (Stoughton et al. 2002), we 



can test whether the scaling relations show similar behav- 
ior for systems with smaller fluxes more typical of galaxy 
groups. A complementary approach (that we do not ap- 
ply here) to understanding scaling relations for low-mass 
systems is to identify the systems optically and use stack- 
ing analysis to measure their ensemble properties. This 
approach has been applied previously with RASS data 
(Dai et al. 2007; Rykoff et al. 2008). 

We describe the data and the group sample in § 2. 
In § 3, we review the caustic technique and use it to 
estimate the group mass profiles, discuss cluster scal- 
ing relations, and compare the caustic mass profiles to 
simple parametric models. We discuss some individ- 
ual groups in §4. We compare the caustic mass pro- 
files to X-ray and virial mass estimators in § 5. We 
discuss our results and conclude in §6. We assume 
H = 100 h km s^/i^Mpc -1 ,^ = 0.3, fl A = 0.7 
throughout. 

2. THE 400D-SDSS GROUP SAMPLE 
2.1. Sloan Digital Sky Survey 

The Sloan Digital Sky Survey (SDSS, Stoughton et al. 
2002) is a wide-area photometric and spectroscopic sur- 
vey at high Galactic latitudes. The Fifth Data Release 
(DR5) of SDSS includes 8000 square degrees of imag- 
ing data and 5740 square degrees of spectroscopic data 
(Adelman-McCarthy et al. 2007). 

The spectroscopic limit of the main galaxy sample of 
SDSS is r=17.77 after correcting for Galactic extinction 
(Strauss et al. 2002). CAIRNS and CIRS found that in- 
fall patterns were detectable in clusters sampled to about 
M* + l, or z <0.1 for SDSS data. For the current sample, 
the X-ray fluxes are much smaller than those of CAIRNS 
or CIRS, so these groups may not be as well sampled. If 
the 400d groups are less massive on average, they should 
contain fewer luminous galaxies than the CAIRNS and 
CIRS clusters. We therefore expect that infall patterns 
may be less common and/or poorly sampled in the 400d 
groups. 

Note that SDSS is ^85-90% complete to the nominal 
spectroscopic limit. The survey has «7% incomplete- 
ness due to fiber collisions (Strauss et al. 2002), which 
are more likely to occur in dense cluster fields. Because 
the target selected in a fiber collision is determined ran- 
domly, this incompleteness can theoretically be corrected 
for in later analysis. From a comparison of SDSS with the 
Millennium Galaxy Catalogue, Cross et al. (2004) con- 
clude that there is an additional incompleteness of ^7% 
due to galaxies misclassified as stars or otherwise missed 
by the SDSS photometric pipeline. For our purposes, 
the incompleteness is not important provided sufficient 
numbers of group galaxies do have spectra. 

2.2. X-ray Surveys for Clusters and Groups 

In CIRS, we studied clusters contained in catalogs 
based on the ROSAT All-Sky Survey (RASS Voges et al. 
1999). RASS is a shallow survey but it is sufficiently deep 
to include nearby, massive clusters and groups. RASS 
covers virtually the entire sky and is thus the most com- 
plete X-ray cluster survey for nearby clusters and groups. 
The flux limits of RASS-based surveys are «3xl0~ 12 crg 
s^cm- 2 (Ebeling et al. 1998; Bohringer et al. 2000a). 

The wide field-of-view of the ROSAT PSPC instru- 
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ment provides a large area for serendipitous discovery 
of sources in pointed observations. Several cluster cat- 
alogs have been created with this goal. The largest of 
these catalogs is the 400d survey, which contains essen- 
tially all pointed observations suitable for extragalactic 
surveys (Burenin et al. 2007), yielding a total survey of 
about 400 square degrees. Extended X-ray sources are 
detected with a wavelet technique and fluxes are mea- 
sured in the ROSAT band (0.1-2.4 kcV). The flux limit 
of the 400d survey is «1.4x 10~ 13 erg s _1 cm -2 , a factor 
of ^20 fainter than the RASS-based catalogs. The ef- 
fective area of the survey depends on the flux limit and 
is calibrated with simulated observations of real ROSAT 
cluster data. The detection algorithm is optimized for 
luminous clusters (Lx > 5 x 10 43 erg s _1 at moderate 
rcdshift (z > 0.3). The groups studied here have smaller 
X-ray luminosities and are thus expected to be physically 
smaller than higher-luminosity clusters, but they are also 
at much lower redshifts, so their angular sizes are not 
too different from the luminous clusters at z>0.3. The 
ROSAT observations are typically not deep enough to 
separate diffuse group emission from contributions from 
individidual galaxy halos. 

We restrict our analysis to clusters and groups with 
z <0.10 (§2.1). Our sample contains 16 groups within 
the SDSS DR5 spectroscopic footprint. We will refer to 
this sample as the 400d-SDSS groups hereafter. Figure 1 
shows the X-ray luminosities of these systems compared 
to CIRS. The range of Lx in the two surveys is similar, 
but the 400d-SDSS sample contains many more systems 
with small Lx than does CIRS. 

3. RESULTS 
3.1. Infall Patterns around X-ray Groups 

We first search for well-defined infall patterns around 
the 400d-SDSS groups. Analogously to CIRS, we plotted 
radius-rcdshift diagrams for all groups in the 400d X-ray 
catalog covered by DR5 with z <0.10. We assign a "by- 
eye" classification of each group's infall pattern: "clean" 
for groups with few background and foreground galaxies, 
"intermediate" for groups with apparent infall patterns 
but significant contamination from either related large- 
scale structure or foreground and background objects, 
and "weak" for groups with little apparent infall pat- 
tern. CIRS successfully applied the caustic technique to 
clusters classified as "intermediate" or "weak" ; the clas- 
sification scheme is thus fairly conservative. We use this 
classification scheme only to show the dependence of the 
infall pattern appearance on group mass (using Lx as a 
proxy) and the sampling depth. 

Figure 1 shows the dependence of this subjective clas- 
sification on X-ray luminosity and redshift. The 400d- 
SDSS sample contains 16 groups; six contain "clean" in- 
fall patterns and three (19%) show weak infall patterns. 
The percentage in the last category is larger than for 
the CIRS clusters (4%), although the significance of this 
difference is difficult to assess (the 400d-SDSS sample is 
small and the categories are not robustly defined). 

Figure 1 demonstrates the expanded parameter space 
covered by the 400d-SDSS groups compared to the 
CIRS sample. In particular, the 400d-SDSS sample in- 
cludes many more systems with Lx^l0 42 /i~ 2 crg s _1 
than CIRS. Table 1 lists the groups in the 400d-SDSS 
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Fig. 1. — Redshift versus X-ray luminosity (0.1-2.4 keV) for 
X-ray clusters from CIRS (small squares) and 400d-SDSS groups 
(large hexagons and stars) contained in the SDSS DR5 spectro- 
scopic survey region. Filled hexagons, open hexagons and six- 
pointed stars indicate 400d-SDSS groups with "clean" , "intermedi- 
ate" , and "weak" infall patterns respectively. These classifications 
for CIRS clusters are represented by filled squares, open squares, 
and crosses. Large open circles around some of the 400d-SDSS 
groups indicate systems contaminated by nearby structure (see 
§3.3). The dashed line and lower solid lines show flux limits of 
fx >1.4 and 3 xl0~ 13 erg s" 1 respectively; the effective survey 
areas at these limits for the 400d survey are 232 and 370 square 
degrees respectively (not all of this area is contained within the 
SDSS DR5 spectroscopic footprint). The upper solid line shows 
the flux and redshift limits of the CIRS cluster sample. The dash- 
dotted line shows the rcdshift and luminosity limits of CAIRNS. 

sample, their X-ray positions and luminosities, their cen- 
tral redshifts and velocity dispersions (see below), and 
the projected radius i? comp within which the SDSS DR5 
spectroscopic survey provides complete spatial coverage. 
For several groups, the caustic pattern disappears be- 
yond Rcomp because of this edge effect. Figures 2-3 show 
the infall patterns for the 400d-SDSS sample. 

Another factor which determines the presence or ab- 
sence of a "clean" infall pattern is the surrounding large- 
scale structure. For example, the redshift-radius dia- 
grams of groups within the infall regions of more massive 
clusters reflect the kinematics of the cluster's infall region 
(e.g., A2199, see Rines et al. 2002). Closer inspection of 
the 400d-SDSS groups reveals that the caustic technique 
misidentifies four of the 16 groups. That is, the galaxy 
group associated with the 400d X-ray source is embedded 
in a larger system containing more galaxies. The hierar- 
chical clustering algorithm used by the caustic technique 
to define the center of the system then mistakes the cen- 
ter of the larger system as the center of the 400d group. 
Ignoring these four groups, the percentage of "weak" in- 
fall patterns increases from 19% to 25% (three of the 12 
remaining 400d-SDSS groups). Although the statistics 
are poor, this result suggests that the 400d-SDSS groups 
might have less well-defined infall patterns than the gen- 
erally more massive CIRS clusters. 

We attempt to quantify this difference more robustly 
by defining a simple statistic to quantify the contrast of 
the infall pattern with respect to the local background. 
We define the contrast C200 to be the ratio of the number 



4 



Rines and Diaferio 



B 

CO 

o 



N 
O 
< 



4 
2 


-2 
4 
2 


-2 
4 
2 


-2 
4 
2 


-2 




TT 



4 

-51 oT" 



6 



*°4 

°°o o 8 



cl0900 - 4£>4» 




h — I — I — M- 



• o" 



C11033 



H — I — — I — 1 8 I "J °i" J-. J 




cll058 ° °»° 8 • o* 8 

J I I J I t J-2 I i ° i 



J «l_ 



2 4 

1 I 1 a. 



6 

-|— IT 



" 1 ' I • » oic 



V 

o o_ ° 



cIIOIO 






■■ — 1—4 — i — 1 







■ f.1 *r- 





— 1 — 1 — 

o 


o o o 








o o 



o- 
o 


'-_ 1" 1 J — 


o 




o 



K4 L ° ■• 



cll039 



H — I — h 



goo e „ S °°oo o o ° 



H — F 



° • • „• *o°° ^o" ° <#!.. „-*-~^ • V Oo"*" * 



o o ° o o -S 



C11159 



^ -2 



- -2 



2 4 60 2 4 6 

Rpth- 1 Mpc) 

Fig. 2. — Rcdshift versus radius for SDSS galaxies around the first eight X-ray groups in the 400d-SDSS sample. The caustic pattern 
is evident as the trumpet-shaped regions with high density. The solid lines indicate our estimate of the location of the caustics in each 
group. Vertical lines in each panel indicate the radius R comp where the spatial coverage of the SDSS DR5 spectroscopic survey is no longer 
complete (no line is shown if Rcomp > 7/i _1 Mpc). Squares around group names indicate systems with contamination from nearby groups 
or clusters (see §3.3). Figure 3 shows similar plots for the rest of the sample. 



of galaxies within projected radius r2oo within the caus- 
tics N mem to the number of "near-background" galaxies 
Nbkgdi those within ±5(T2oo of the cluster rcdshift but 
outside the caustics. The radius r 2 oo is determined from 
the caustic mass profile (see below) and 0200 is the pro- 
jected velocity dispersion of galaxies within the caustics 
and within the projected radius r 2 oo- Using this defini- 
tion of contrast, there is no significant correlation of C200 
with cluster mass for the CIRS clusters. Further, there 
is no significant difference between the distributions of 
the C200 statistic between the 400d-SDSS (excluding the 
four problematic groups) and CIRS samples (the sam- 
ples differ at 59% confidence using a K-S test). This 
result suggests that any difference between the contrasts 
of the two samples is too subtle for our simple contrast 
statistic. We caution the reader that this contrast statis- 
tic is not well tested and may have subtle dependencies 
on the caustic method used to define membership. It 



does provide a simple, quantitative way of comparing 
the 400d-SDSS sample to the CIRS sample. 
We discuss the individual groups in more detail in §3.3. 

3.2. Caustics and Mass Profiles 

CAIRNS and CIRS showed that caustic masses of clus- 
ters agree well with mass estimates from both X-ray 
observations and Jeans' analysis at small radii, where 
dynamical equilibrium is expected to hold (Rines et al. 
2003, CIRS). Lokas & Mamon (2003) confirm that the 
mass of Coma estimated from higher moments of the 
velocity distribution agrees well with the caustic mass 
estimate (Geller et al. 1999). Recently, Diaferio et al. 
(2005) showed that caustic masses agree with weak lens- 
ing masses in three clusters at moderate redshift. 

We briefly review the method DG and D99 developed 
to estimate the mass profile of a galaxy cluster by iden- 
tifying caustics in redshift space. The method assumes 
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Fig. 3.— See Figure 2. 



that clusters form in a hierarchical process. Applica- 
tion of the method requires only galaxy redshifts and sky 
coordinates. Toy models of simple spherical infall onto 
clusters produce sharp enhancements in the phase space 
density around the system. These enhancements, known 
as caustics, appear as a trumpet shape in scatter plots 
of redshift versus projected clustercentric radius (Kaiser 
1987). DG and D99 show that random motions smooth 
out the sharp pattern expected from simple spherical in- 
fall into a dense envelope in the redshift-projected radius 
diagram (see also Vedel & Hartwick 1998). The edges 
of this envelope can be interpreted as the escape veloc- 
ity as a function of radius. Galaxies outside the caustics 
are also outside the turnaround radius. The caustic tech- 
nique provides a well-defined boundary between the infall 
region and interlopers; one may think of the technique as 
a method for defining membership that gives the cluster 
mass profile as a byproduct. In fact, the caustic tech- 
nique can be used as a membership classification for any 
gravitationally bound system. For instance, Serra et al. 



(2009) used caustics to identify interloper stars in five 
dwarf spheroidal galaxies, and Brown et al. (2009) use 
caustics to identify stars not bound to the Milky Way. 

Operationally, we identify the caustics as curves which 
delineate a significant decrease in the phase space density 
of galaxies in the projected radius-redshift diagram. The 
details of the caustic technique used here are identical to 
those described in CIRS, with the following exception. 
Because 400d-SDSS groups are expected to be less mas- 
sive than CIRS clusters, we isolate the groups initially 
by studying only galaxies within R p <7/i _1 Mpc and 
±4000km s _1 of the nominal group centers from the X- 
ray catalogs (compared to 10/i _1 Mpc and ±5000km s _1 
for the CIRS clusters). We perform a hierarchical struc- 
ture analysis to locate the centroid of the largest system 
in each volume. This analysis sometimes finds the center 
of another system in the field. In these cases, limiting the 
galaxies to a smaller radial and/or redshift range enables 
the algorithm to center on the desired group. For some 
groups, no satisfactory match between the hierarchical 



6 



Rines and Diaferio 



TABLE 1 
400D-SDSS Basic Properties 
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Note. — Col. (1) lists the number of the group in the 400d catalog (Burenin et al. 2007). Col. (2) lists the name of the group following 
the convention of Burenin et al. (2007). Col. (8) gives our subjective classification of the contrast between the group and foreground and 
background structure. A Flag of 2 indicates a "clean" infall pattern with little contaminating structure, while a Flag of indicates a 
"weak" infall pattern for groups where the infall pattern is weakly present and/or significantly contaminated by nearby structure (see text 
for details). Col. (9) gives the projected radius of the nearest edge of the footprint of the SDSS spectroscopic survey. 

a Strongly affected by a nearby group or cluster. See §3.3 for 
details. 

b Burenin et al. (2007) list the rcdshift as 2=0.07. 



center and X-ray center is possible; for these groups we 
impose the X-ray center on the analysis (Table 2). 

We discuss the groups individually in §3.3. Figures 2-3 
show the caustics and Figures 4-5 show the associated 
mass profiles. Note that the caustics extend to different 
radii for different groups. D99 shows that the appear- 
ance of the caustics depends strongly on the line of sight; 
projection effects can therefore account for most of the 
differences in profile shape in Figures 2-3 without invok- 
ing non-homology among clusters. We use the caustics 
to determine group membership. Here, the term "group 
member" refers to galaxies both in the virial region and 
in the infall region and inside the caustics. Figures 2-3 
show that the caustics effectively separate group mem- 
bers from background and foreground galaxies, although 
some interlopers may lie within the caustics. Note that, 
for some groups, the small samples of galaxies yield unre- 
alistically low estimates of the uncertainties in the caustic 
mass profiles. 

The D99 algorithm we use to identify the caustics gen- 
erally agrees with the lines one would draw based on a 
visual impression. This consistency suggests that sys- 
tematic uncertainties in the caustic technique are domi- 
nated by projection effects rather than the details of the 
algorithm. We now discuss the individual groups in more 
detail. 

3.3. Comments on Individual Groups 

Groups share many common features, but each system 
is unique. We describe the most relevant aspects of each 
400d-SDSS group below. 

• cl0810+4216 The X-ray emission is centered on a 
r=14.2 (M r «-22.3) galaxy. The hierarchical cen- 



ter is located at the position of a strongly bound 
galaxy pair 387/i _1 kpc to the NE. The X-ray center 
corresponds to the position of the Brightest Group 
Galaxy (BGG) and is probably the dynamical cen- 
ter of this system. The infall diagram produced 
when adopting the X-ray center as the system cen- 
ter suggests that four sub-L* galaxies have slightly 
higher redshifts than the bulk of the group mem- 
bers. Excluding these galaxies from the caustic- 
selected sample would reduce the inferred velocity 
dispersion by ~20%. The caustic technique classi- 
fies these four galaxies as members using either the 
X-ray or the hierarchical center. 

• cW820+5645 The X-ray emission is centered on 
a pair of interacting galaxies with r=15.7 and 
r=f6.0. A plot of group members on the sky re- 
veals two main clumps of galaxies, one on the X- 
ray center and the other approximately centered 
on a r = 14.6 galaxy (MCG +09-14-024) without 
a redshift. This system is possibly a group-group 
merger. 

• cl0900+3920 The X-ray emission is centered <V 
from a pair of galaxies with r=15.9 and r=17.1. 
The Brightest Group Galaxy (r=15.4) lies ~4' 
north of the X-ray center. 

• cllOlO+5430 The X-ray emission is centered on a 
r=14.2 galaxy, UGC 06057. This galaxy has no 
SDSS rcdshift, but a redshift from Carballo et al. 
(1995) places it in the group. 

• cll033+5703 The X-ray emission is centered on a 
r=14.2 galaxy. The spatial distribution of galax- 
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Fig. 4. — Comparison of caustic mass profiles to those estimates from the virial theorem and the projected mass estimator. The thick 
solid lines show the caustic mass profiles and the thin lines show the lcr uncertainties in the mass profiles (for some systems, the thin lines 
blend with the thick lines). The axes are identical in all panels. The vertical bars indicate r2oo an d the maximum radius of the caustic mass 
profile (the smaller of r max , the extent of the infall pattern, and rt, the turnaround radius). Vertical dashed lines indicate rt for groups 
where the infall pattern truncates before rt- Red and green shaded regions show the formal lcr uncertainties in the virial and projected 
mass profiles. Squares around group names indicate groups with contamination from nearby groups or clusters (see §3.3). Figure 5 shows 
similar plots for the rest of the sample. 



ies near this group shows that the X-ray group 
is located ^800/i _1 kpc from a group of galaxies 
at the same redshift centered on CGCG 290-048 
(r=13.5), a radio galaxy with a jet 1.3/i _1 Mpc in 
extent (Masson 1979). The 400d ROSAT point- 
ing does not include the CGCG 290-048 group, 
but it has been classified as an optically-selected 
group by Miller et al. (2005) and Merchan & Zan- 
divarez (2005). The dynamics of cll033+5703 are 
probably perturbed by the proximity of the CGCG 
290-048 group (for an analogous configuration, see 
A2197E/W and A2199 Rines et al. 2002). Indeed, 
the CGCG 290-048 group displays a more symmet- 
ric infall pattern than does cll033+5703, suggest- 
ing that the former system dominates the dynamics 
of the two systems. 



• cll039+3947 The X-ray emission is approximately 
centered on a r=15.0 galaxy. A bright (r=15.2) 
spiral galaxy with several fainter companions lies 
about 2/i~ 1 Mpc NW of the group center. This lat- 
ter system may be an infalling group. 

• cll058+0136 The X-ray emission is centered on 
UGC 06057, a r=13.5 galaxy without a SDSS red- 
shift but confirmed as a group member using a 
redshift from Smith et al. (2000). This group is 
6.'6 from the optical position of Abell 1139 (Abell 
et al. 1989) and previously has been identified with 
A1139 by Ebeling et al. (1998) and Bohringer et al. 
(2000b). 

• clll59+5531 The X-ray emission is centered on a 
r=14.1 (M r Ps-23.0) galaxy. This group was clas- 
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Fig. 5. — See Figure 4. 



sified as an X-ray overluminous elliptical galaxy 
(OLEG) by Vikhlinin et al. (1999). 

cll227+0858 The X-ray emission is centered on a 
r=15.0 galaxy without a SDSS rcdshift. A rcdshift 
from Slinglend et al. (1998) shows this galaxy to be 
a member. The group lies ~0.7/i Mpc from the 
Brightest Cluster Galaxy (BCG) of Abell 1541, a 
cluster at the same redshift. The hierarchical cen- 
ter is coincident with the core of A 1541, so we in- 
stead use the X-ray center to estimate the caustics. 
Bohringer et al. (2000b) report an X-ray flux of 
(2.5±0.5)xl0- 12 erg s^cm -2 for A1541, ~6 times 
larger than the flux from cll227+0858. This differ- 
ence suggests that A 1541 dominates the dynamics 
of cll227+0858 (see cll033+5703 above). The dy- 
namical parameters reported for cll227+0858 thus 
partially reflect the properties of A1541, a more 
massive and more X-ray luminous system. Consis- 
tent with this idea, 7 galaxies within 0.4/i -1 Mpc 
of the X-ray center have a velocity dispersion 



(594± 2l J 2 km s x ) 20% smaller than the value re- 
ported in Table 1; at slightly larger radii there is a 
large spike in redshift space due to A1541 (Figure 
3). 



• cll236+1240 This extended X-ray source is ap- 
proximately centered on IC 3574, a r=14.1 (M r «- 
22.4) galaxy at 0=0.067 (Burenin et al. 2007, no 
SDSS spectrum available). SDSS spectra, however, 
show that there are 15 galaxies at zwO.044 within 
0.5/i _1 Mpc of the X-ray center, while there are no 
other galaxies within this radius at zwO.067. The 
brightest of these galaxies is r=15.5 (Af r ss-20.2). 
The X-ray emission is therefore associated with a 
single bright galaxy (similar to a fossil group, al- 
though most fossil groups do contain faint galaxies 
at the same redshift) or with a group of less lu- 
minous galaxies (or some combination). We adopt 
the latter interpretation that cll236+1240 is asso- 
ciated with the group of less luminous galaxies. 
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TABLE 2 

400D-SDSS Hierarchical Centers and Offsets From X-ray Centers 



Croup 


Hierarchical Center 
RA (J2000) DEC (J2000) 


Acz" 1 
km s 


AR & 

It KpC 


ciuoiu-r^t^io 


122.74532 


42.33660 


Ulo 


OO f 


ClUoZU-rDO^iO 


125.11101 


56.75749 






^•inonn_i_Qoon 
ciuyuu-roy^u 


135.02003 


39.34003 






CI 1U IUtJ^OU 


152.55563 


54.55338 


-ouu 


11^ 
11J 


cilUoo-t-o 1 Uo 


158.46603 


57.05308 






cll039+3947 


159.87998 


39.79392 






C11058+0136 


164.53830 


01.62267 


-393 


30 


ell 159+5531 


179.96301 


55.53222 






cll227+0858 b 


186.87500 


08.82286 


-34 c 


_d 


C11236+1240 


189.06404 


12.44854 


002 


149 


cll329+1143 


202.32590 


11.65740 


-17 


89 


cll343+5546 b 


205.90356 


55.62635 


-436 


518 


cll533+3108 


233.32102 


31.14862 






cll630+2434 


247.59420 


24.57920 


702 


92 


cll631+2121 


247.86281 


21.51073 


831 


_d 


cl2137+0026 


324.25590 


00.40342 


127 


121 



a No value listed for the 6 groups where the X-ray center is used 
for the caustic analysis. 

b Strongly affected by a nearby group or cluster. See §3.3 for 
details. 

c Burcnin ct al. (2007) list a rcdshift of z=0.087 and reference 
Burke et al. (2003). However, Burke ct al. (2003) list a redshift of 
2=0.090 for this group. 

d X -ray center adopted for spatial position and hierarchical 
central redshift. The hierarchical centers are offset on the 
sky by Ai?=685 /i^kpc and 753/i" 1 kpc for cll227+0858 and 
C11631+2121 respectively. 

• cll329+1143 The X-ray emission is centered 
between NGC 5179 (r=13.3) and NGC 5171 
(r=12.9). The latter has no SDSS spectrum, but a 
redshift from Falco et al. (1999) places NGC 5171 
in the group. This group is also known as MKW 
11 and its infall region was previously analyzed in 
CIRS. Note that the X-ray flux from the 400d cat- 
alog [(1.3±0.2)xl0~ 12 crg s _1 cm -2 ] is significantly 
smaller than the flux from RASS data Bdhringer 
et al. (2000b) used in CIRS [(5.7±0.8) x 10~ 12 erg 
s _1 cm -2 ]. The smaller flux from the 400d catalog 
would place this system below the flux limit used 
in CIRS. Osmond et al. (2004) analyze an XMM- 
Newton observation of this group and determine 
a flux of (2.5±0.4)xl0~ 12 erg s^cm" 2 in the (0.5- 
2.0 kcV) band, intermediate between the RASS and 
400d values. A possible explanation of the varying 
flux estimates is that a background group south- 
east of NGC 5171 (centered on a red galaxy with 
SDSS photometric redshift z=0.21±0.01) may have 
been included in the RASS flux, while the 400d 
flux may exclude some emission associated with 
group members. Osmond et al. (2004) further find 
evidence for interactions between the group mem- 
bers NGC 5171 and NGC 5176 and the surround- 
ing gas. We find a slightly smaller velocity disper- 
sion (363^33 km s _1 ) than do Osmond & Ponman 
(2004, 494±99 km s _1 ), who noted that their value 
of a p was larger than expected based on the mea- 
sured Tx=0.96±0.04 keV. Our smaller value of a p 
mitigates this discrepancy. The NGC 5129 X-ray 
group lies ^3h~ 1 Mpc northeast of cll329+1143 at 
a similar redshift; it is possible that these groups 
are physically associated. 



• cll343+5546 The X-ray emission is centered 
on a r=15.8 galaxy. The group is located 
llf4«500/i" 1 kpc from the optical center of A1783. 
The hierarchical center is coincident with the core 
of A1783, so the reported dynamical properties 
reflect those of the Abell cluster and not the 
400d group. Unlike cll227+0858/Al541 above, 
cll343+5546 and A1783 are too close together 
to obtain a mass profile for cll343+5546 alone. 
Ledlow et al. (2003) report an X-ray flux of 
(0.51±0.12)xl0- 12 crg s^cm- 2 for A1783, -2.5 
times larger than the flux of cll343+5546. The 
galaxies in A1783 exhibit significant substructure, 
perhaps indicating that the cluster is not yet re- 
laxed. 

• cll533+3108 The X-ray emission is centered on a 
r=15.3 galaxy possibly undergoing a merger with 
a r=16.7 galaxy; neither galaxy has a SDSS spec- 
trum. A redshift from Small et al. (1997) places 
the brighter galaxy in the group. The central 
galaxy pair is surrounded by several bright galax- 
ies. This group is also known as Abell 2092. Wu 
et al. (1999) list a velocity dispersion of a v — 
504lgg 5 kms -1 , compared with our estimate of 

(T p = 421^51 km s _1 . Our estimate is based on 
a larger number of galaxies. 

• cll630+2434 The X-ray emission is centered be- 
tween a r=14.6 and a r=15.2 galaxy. The brighter 
galaxy (IC 4607) has no redshift available from ei- 
ther SDSS or the literature. 

• cll631+2121 The X-ray center lies among seven 
red galaxies with r=16. 5-18.0 (one a spectroscop- 
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ically confirmed member); a r—l6A spiral galaxy 
(and spectroscopically confirmed member) lies ~1' 
northeast of the X-ray center. The hierarchical cen- 
ter is ~750/i -1 kpc northeast of the X-ray center, 
and <~1.'5 southwest of a r=15.4 group member. 
Many group members lie close to the line connect- 
ing the X-ray and hierarchical centers. We center 
the mass profile on the X-ray center and the hierar- 
chical redshift (Table 2). Additional redshifts, es- 
pecially of the galaxies near the X-ray center, would 
help clarify the nature of this system. 

• cl2137+0026 The X-ray emission is centered be- 
tween a r=14.3 early- type galaxy and a r=14.4 
spiral galaxy. Both are group members, although 
the early-type galaxy has no SDSS spectrum (Falco 
et al. 1999). 

3.4. Virial and Turnaround Masses and Radii 

The caustic mass profiles allow direct estimates of the 
virial and turnaround radius in each group. For the virial 
radius, we estimate r 2 oo i r A is the radius within which 
the enclosed average mass density is Ap c , where p c is the 
critical density) by computing the enclosed density pro- 
file [p(< r) = 3M(< r)/47rr 3 ]; r 2 oo is the radius which 
satisfies p(< r 2 oo) = 200 p c . In our adopted cosmology, 
a system should be virialized inside the slightly larger 
radius ^rioo ~ l-3r 2 oo (Eke et al. 1996). We use r 2 oo be- 
cause it is more commonly used in the literature and thus 
allows easier comparison of results. For the turnaround 
radius r t , we use equation (8) of Regos & Geller (1989) 
assuming Q m — 0.3. For this value of Q TO , the enclosed 
density is 3.5p c at the turnaround radius. Varying TO 
in the range 0.02-1 only changes the inferred value of r t 
by ±10%; the uncertainties in r t from the uncertainties 
in the mass profile are comparable or larger (D99; Rines 
et al. 2002). If the w parameter in the equation of state 
of the dark energy (Pa = wp^) satisfies w > — 1, the 
dark energy has little effect on the turnaround overden- 
sity (Gramann & Suhhoncnko 2002). 

Table 3 lists r 20 o, r tl and the masses M 20 o and M t en- 
closed within these radii. For cl0820+0945, cll329+1143, 
and cll631+21214, the maximum extent of the caustics 
fmax is smaller than r t . For these groups, r t is a min- 
imum value assuming that there is no additional mass 
outside r max . The best estimate of the mass contained 
in infall regions clearly comes from those groups for which 
r m ax > Tt- The average mass within the turnaround ra- 
dius for these groups is 2.4±0.4 times the virial mass 
M200, demonstrating that groups are still forming in 
the present epoch and consistent with the CIRS esti- 
mate of 2.19±0.18. Similarly, the average turnaround ra- 
dius is (5.01±0.25)r 2 oo for groups with r max > r t , again 
in agreement with the CIRS result of (4.96±0.08)r 20 o- 
Simulations of the future growth of large-scale struc- 
ture (Gramann & Suhhonenko 2002; Nagamine & Loeb 
2003; Busha et al. 2003) for our assumed cosmology 
(tt m — 0.3, — 0.7) suggest that galaxies currently in- 
side the turnaround radius of a system will continue to be 
bound to that system. Our results for the turnaround ra- 
dius and mass agree with the predictions of Busha et al. 
(2005), who find that the ultimate mass of dark mat- 
ter haloes in simulations is 1.9M 2 oo- The agreement be- 
tween the estimates of M t /M 20 o and r t /V 2 oo between the 



CIRS and 400d-SDSS samples suggests that the overall 
shapes of cluster mass profiles into their infall regions are 
not strongly dependent on cluster mass into the group 
regime. We compare the virial masses M 20 o from the 
caustics with masses calculated using the virial theorem 
in §3.5. 

One striking result of this analysis is that the caus- 
tic pattern is often visible beyond the turnaround radius 
of a group (similar to some CIRS clusters). This re- 
sult suggests that even groups may have strong dynamic 
effects on surrounding large-scale structure beyond the 
turnaround radius. For our assumed cosmology, this 
large-scale structure is probably not bound to the group. 

3.5. Comparison to Virial and Projected Mass 
Estimates 

We apply the virial mass and projected mass estima- 
tors (Heisler et al. 1985) to the 400d-SDSS groups. For 
the latter, we assume the galaxies are on isotropic orbits. 
We must define a radius of virialization within which the 
galaxies are relaxed. We use r 20 o (Table 3) and include 
only galaxies within the caustics. We thus assume that 
the caustics provide a good division between group galax- 
ies and interlopers (see Figures 2-3). For details on these 
mass calculations, see CIRS. Table 4 lists the virial and 
projected mass estimates. 

Figure 6 compares the virial and caustic mass esti- 
mates at r 2 oo- The mean ratios of these estimates are 
M v /M c — 0.91 ± 0.12. The caustic mass estimates are 
consistent with virial mass estimates even assuming a 
correction factor C ~ 0.1 — 0.2M vir for the surface pres- 
sure, consistent with the best-fit NFW profiles (see also 
Carlberg et al. 1997; Girardi et al. 1998; Koranyi & Geller 
2000; Rines et al. 2003, CIRS). Note that Lemze et al. 
(2009) suggest a possible modification of the form factor 
used in the caustic technique based on a comparison 
of the caustic and lensing mass profiles of A1689. Lopes 
et al. (2009) find a smaller scatter in scaling relations us- 
ing virial masses rather than caustic masses for the CIRS 
clusters. 

Figures 4-5 compare the mass profiles estimated from 
the caustics, virial theorem, and projected mass estima- 
tor. The virial and projected mass profiles M{< Rp) are 
calculated using all galaxies within the projected radius 
Rp (note that the virial assumptions hold only at the 
virial radius) . The projected mass estimator consistently 
overestimates the mass at small radii and underestimates 
the mass at large radii relative to the other profiles. This 
behavior suggests that this estimator is best for estimat- 
ing virial masses but not mass profiles. The virial and 
caustic mass profiles generally agree although there are 
many groups with large disagreements. The caustic mass 
profiles do not appear to consistently overestimate or un- 
derestimate the mass relative to the virial mass profiles. 
This result supports our use of caustic mass profiles as a 
tracer of the total group mass profile in the next section. 

3.6. The Shapes of Group Mass Profiles 

We fit the mass profiles of the 400d-SDSS groups to 
three simple analytic models. The simplest model of 
a self-gravitating system is a singular isothermal sphere 
(SIS). The mass of the SIS increases linearly with radius. 
Navarro et al. (1997) and Hernquist (1990) propose two- 
parameter models based on CDM simulations of haloes. 
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TABLE 3 

400D-SDSS Characteristic Radii and Masses 



Group 


f500 


»*200 


n 


Tmax 


^200 


M t 


M t /M 20 o 






fc _1 Mpc 


h _1 Mpc 


h _1 Mpc 


1O 14 M 


1O 14 M 




cl0810+4216 


0.52 


0.83 


4.05 


9.49 


1.32± 


0.05 


2.70±0.14 


2.05 


cl0820+5645 a 


0.67 


1.09 


>4.55 


3.94 


3.03± 


0.02 


3.84±0.02 


1.27 


cl0900+3920 


0.35 


0.58 


4.17 


9.19 


0.45± 


0.06 


2.95±0.48 


6.58 


cllOlO+5430 


0.20 


0.42 


2.66 


9.80 


0.17± 


0.04 


0.77±0.25 


4.45 


cll033+5703 a 


0.38 


0.60 


3.14 


9.60 


0.51± 


0.02 


1.26±0.06 


2.48 


cll039+3947 


0.29 


0.45 


2.08 


3.64 


0.21± 


0.04 


0.37±0.08 


1.74 


cll058+0136 


0.51 


0.86 


4.17 


>10 


1.51± 


0.00 


2.96±0.00 


1.96 


clll59+5531 


0.39 


0.63 


2.86 


9.29 


0.57± 


0.00 


0.95±0.00 


1.67 


cll227+0858 a 


0.56 


0.97 


4.09 


4.34 


2.13± 


0.83 


2.79±1.16 


1.31 


cll236+1240 


0.49 


0.71 


2.88 


9.70 


0.84± 


0.00 


0.97±0.00 


1.16 


cll329+1143 


0.40 


0.59 


>2.66 


1.82 


0.48± 


0.00 


0.77±0.00 


1.58 


cll343+5546 a 


0.53 


0.90 


4.70 


>10 


1.69± 


0.06 


4.23±0.20 


2.50 


cll533+3108 


0.32 


0.54 


2.50 


2.63 


0.36± 


0.23 


0.64±0.47 


1.76 


cll630+2434 


0.52 


0.76 


3.12 


9.09 


1.03± 


0.18 


1.24±0.23 


1.20 


cll631+2121 


0.40 


0.60 


>2.72 


1.82 


0.49± 


0.39 


0.82±0.71 


1.67 


cl2137+0026 


0.35 


0.51 


2.67 


8.99 


0.32± 


0.01 


0.78±0.05 


2.45 



Note. — Cols. (l)-(3) list r$ for three values of ovcrdcnsity <5, where r$ is the radius within which the enclosed mass density is <5 times 
the critical density (we adopt <5=3.5 as the turnaround ovcrdcnsity). Col. (4) gives the turnaround radius r t for each group. Col. (5) gives 
the maximum radius r max within which the caustics arc detected. Note that when r max < rt, the estimate of rt is an underestimate if 
any mass lies between r ma x and rt. 

a Strongly affected by a nearby group or cluster. See §3.3 for 
details. 



TABLE 4 

400D-SDSS Virial and Projected Masses 



Group 


r200 


M'200 








fe-'Mpc 


1O 14 M 


10 14 M Q 


1O 14 M 


cl0810+4216 


0.83 


1.32± 


0.05 


0.79±0.19 


1.14±0.17 


cl0820+5645 a 


1.09 


3.03± 


0.02 


3.15±0.65 


3.15±0.41 


cl0900+3920 


0.58 


0.45± 


0.06 


0.21±0.06 


0.27±0.04 


cllOlO+5430 


0.42 


0.17± 


0.04 


0.20±0.06 


0.12±0.02 


cll033+5703 a 


0.60 


0.51± 


0.02 


0.45±0.14 


0.35±0.06 


cll039+3947 


0.45 


0.21± 


0.04 


0.08±0.04 


0.18±0.00 


cll058+0136 


0.86 


1.51± 


0.00 


1.02±0.18 


0.76±0.09 


clll59+5531 


0.63 


0.57± 


0.00 


0.24±0.08 


0.41±0.08 


cll227+0858 a 


0.97 


2.13± 


0.83 


1.09±0.21 


3.58±0.45 


cll236+1240 


0.71 


0.84± 


0.00 


0.39±0.16 


0.59±0.12 


cll329+1143 


0.59 


0.48± 


0.00 


0.58±0.12 


0.52±0.07 


cll343+5546 a 


0.90 


1.69± 


0.06 


0.79±0.15 


0.85±0.11 


cll533+3108 


0.54 


0.36± 


0.23 


0.46±0.13 


0.63±0.10 


cll630+2434 


0.76 


1.03± 


0.18 


0.35±0.09 


0.64±0.10 


cll631+2121 


0.60 


0.49± 


0.39 


0.31±0.13 


0.97±0.20 


cl2137+0026 


0.51 


0.32± 


0.01 


0.07±0.03 


0.09±0.02 



Note. — Cols. (4) and (5) give the projected and virial mass estimates for each group computed with all galaxies inside the radius r200 
determined from the caustic mass profile. 

a Strongly affected by a nearby group or cluster. See §3.3 for 
details. 

We note that the caustic mass profiles mostly sample 
large radii and are therefore not very sensitive to the in- 
ner slope of the mass profile. Thus, we do not consider 
alternative models which differ only in the inner slope 
of the density profile (e.g., Moore ct al. 1999). At large 
radii, the best constraints on cluster mass profiles come 
from galaxy dynamics and weak lensing. The caustic 
mass profiles of Coma (Geller et al. 1999), A576 (Rincs 
et al. 2000), A2199 (Rines et al. 2002) and the rest of the 
CAIRNS and CIRS clusters (Rines et al. 2003; Rines & 
Diaferio 2006) provided strong evidence against a singu- 
lar isothermal sphere (SIS) profile and in favor of steeper 
mass density profiles predicted by Navarro et al. (1997, 
NFW) and Hernquist (1990). Only recently have weak 
lensing mass estimates been able to distinguish between 



SIS and NFW density profiles at large radii (Clowe & 
Schneider 2001; Kneib et al. 2003; Mandelbaum ct al. 
2006). 

At large radii, the NFW mass profile increases as ln(r) 
and the mass of the Hernquist model converges. The 
NFW mass profile is 

where a is the scale radius and M(a) is the mass within 
a. We fit the parameter M(a) rather than the character- 
istic density S c (M(a) — 47T(5 c /9 c a 3 [ln(2) — ^] where p c is 
the critical density) because M(a) and a are much less 
correlated than S c and a (Mahdavi et al. 1999). The 
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Fig. 6. — Caustic masses at T200 compared to virial masses at 
the same radius. Symbol types are the same as for Figure 1. Four 
groups that arc strongly affected by nearby systems arc highlighted 
with large circles. Errorbars show la uncertainties and the solid 
line shows a ratio of unity. 

Hcrnquist mass profile is 

M(< r) - M f* (2) 
(r + a H ) 2 

where an is the scale radius and M is the total mass. 
Note that M(a H ) = Ml A. The SIS mass profile is M(< 
r) oc r. We minimize \ and list the best-fit parameters 
a, ^200j the concentration c 2 oo=f"20o/ a : and M 2 oo for the 
best-fit NFW model and indicate the best-fit profile type 
in Table 5. We also list the parameter Cioi=rioi/a; some 
authors prefer to use rioi as the virial radius for this 
cosmology. We perform the fits on all data points within 
the maximum radial extent of the caustics r max listed in 
Table 3 and with caustic amplitude A(r) > 100 km s _1 . 

Because the individual points in the mass profile are 
not independent, the absolute values of % 2 are indicative 
only, but it is clear that the NFW and Hernquist profiles 
provide acceptable fits to the caustic mass profiles; the 
SIS is never the best-fit profile and it is actually excluded 
for nearly all groups. The NFW profile provides a better 
fit to the data than the Hernquist profile for 5 of the 16 
400d-SDSS groups; the remaining 11 are better fit by a 
Hernquist profile. 

Figure 7 shows the shapes of the caustic mass pro- 
files scaled by r 2 oo and M 20 o along with SIS, NFW, and 
Hernquist model profiles. The colored lines show differ- 
ent model mass profiles. The straight dashed line is the 
SIS, the solid lines are NFW profiles with c 20 o=3,5, and 
10, and the curved dashed lines are Hernquist profiles 
with two different scale radii. The average value of c 2u o 
is 6.5, consistent with CIRS and with the values expected 
from simulations for clusters (NFW; Bullock et al. 2001). 
All three model profiles agree fairly well with the caus- 
tic mass profiles in the range (0.1-l)r 2 oo- The SIS only 
fails beyond ~1.5r 2u o; this is why lensing has had trou- 
ble distinguishing between SIS and NFW profiles. As 
discussed in D99, the caustic technique can be subject 
to large variations for individual groups due to projec- 
tion effects. The best constraints on the shapes of cluster 



and group mass profiles are obtained by averaging over 
many lines of sight, or for real observations, over many 
different clusters and groups. 

The concentration parameters c 20 o = r 20 o/a for the 
NFW models are in the range 1-15 (Table 5), in good 
agreement with the predictions of numerical simulations 
(Navarro et al. 1997). Numerical simulations suggest 
that c should decrease slightly (~30%) with increasing 
mass over this mass range, but this decrease is compara- 
ble to the scatter in c in simulated clusters (Navarro et al. 
1997; Bullock et al. 2001). The errorbars in Figure 8 in- 
dicate the average values and la scatter in cioi = rioi/a 
in simulations (Bullock et al. 2001). 

Observations of mass profiles using X-ray mass esti- 
mates (Buote et al. 2007) and weak lensing (Comerford 
& Natarajan 2007) appear to confirm a decrease in c with 
increasing mass, although some massive clusters appar- 
ently have high concentrations (e.g. Lemze et al. 2009). 
In contrast, the CIRS clusters showed a weak positive 
correlation of c with mass (Figure 8). The 400d-SDSS 
groups do not confirm this correlation (Figure 8). A 
Spearman rank-sum test indicates that there is no signifi- 
cant correlation between cioi and Mioi (20% confidence) 
for the 400d-SDSS groups. A Spearman test for the CIRS 
clusters suggests a correlation at the 92.7% confidence 
level. By combining the CIRS and 400d-SDSS data, a 
Spearman test indicates a correlation between concen- 
tration and mass at the 93.6% confidence level (correla- 
tion coefficient 0.19). This result indicates possible ten- 
sion between caustic mass profiles and simulations, as 
the caustic mass profiles indicate that concentrations in- 
crease rather than decrease with increasing cluster mass. 
More data for groups at lower masses would help clarify 
the situation. 

To further explore the dependence of mass profile 
shapes on mass, the right panel of Figure 7 shows the 
scaled mass profiles of the most massive quartile of CIRS 
clusters. Comparing the two panels reveals no obvious 
trend with increasing mass: the profiles of 400d-SDSS 
groups have larger scatter, but the average of the profiles 
appears similar to that of the massive CIRS clusters. 

4. SCALING RELATIONS OF CLUSTERS AND 
GROUPS 

Scaling relations between simple cluster/group observ- 
ables and masses provide insight into the nature of clus- 
ter/group assembly and the properties of various clus- 
ter/group components. Establishing these relations for 
local clusters and groups is critical for future studies of 
clusters in the distant universe with the goal of constrain- 
ing dark energy (Majumdar & Mohr 2004; Lin et al. 
2004). 

We apply the prescription of Danese et al. (1980) to 
determine the mean redshift cz Q and projected velocity 
dispersion a p of each group from all galaxies within the 
caustics. We calculate a p using only the group members 
projected within r 2 oo estimated from the caustic mass 
profile. Note that our estimates of r 2 oo do not depend 
directly on a p . 

One of the simplest X-ray observables of clusters and 
groups is X-ray luminosity. For optical data, the pro- 
jected velocity dispersion a p is a straightforward observ- 
able. In particular, Lopes et al. (2009) reanalyze the 
CIRS clusters and find that the velocity dispersions are 
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TABLE 5 

400D-SDSS Mass Profile Fit Parameters 



Group 


a NFW 


»"200 


C200 


M 20 o 


Best-tit 


ClOl 




h~ 1 Mpc 


h _1 Mpc 




10 14 M© 


Profile 




C10810+4216 


0.185 


0.79 


4.26 


1.14 


N 


5.75 


cl0820+5645 a 


0.205 


0.95 


4.63 


1.98 


H 


6.23 


cl0900+3920 


0.590 


0.60 


1.02 


0.50 


N 


1.49 


cllOlO+5430 


0.415 


0.40 


0.96 


0.15 


H 


1.42 


cll033+5703 a 


0.160 


0.61 


3.83 


0.53 


H 


5.19 


cll039+3947 


0.072 


0.41 


5.75 


0.16 


H 


7.68 


cll058+0136 


0.297 


0.77 


2.59 


1.06 


11 


3.58 


clll59+5531 


0.160 


0.58 


3.64 


0.46 


11 


4.95 


cll227+0858 a 


0.415 


0.75 


1.80 


1.01 


H 


2.59 


cll236+1240 


0.081 


0.62 


7.70 


0.57 


11 


10.20 


cll329+1143 


0.064 


0.60 


9.32 


0.49 


N 


12.30 


cll343+5546 a 


0.200 


0.86 


4.29 


1.46 


N 


5.78 


cll533+3108 


0.138 


0.48 


3.46 


0.25 


H 


4.71 


cll630+2434 


0.051 


0.68 


13.40 


0.75 


11 


17.56 


cll631+2121 


0.092 


0.58 


6.32 


0.47 


H 


8.52 


cl2137+0026 


0.081 


0.50 


6.20 


0.30 


N 


8.27 



Note. — Cols. (2) and (3) give the scale radius ajvFW an <i r 200 of the best-fit NFW profile. Cols. (4) and (6) give the NFW concentration 
parameters eg = rg/apfpyy where rg is the radius within which the enclosed density is 5 times the critical density. 

a Strongly affected by a nearby group or cluster. See §3.3 for 
details. 
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Fig. 7. — Scaled caustic mass profiles for the 400d-SDSS groups (excluding those contaminated by nearby large-scale structure) compared 
to simple models. The thin solid lines show the caustic mass profiles normalized by r2oo and M2oo- The long-dashed line shows a singular 
isothermal sphere, the colored solid lines show NFW profiles (with concentrations C200=3,5,10 from top to bottom at large radii). The 
short-dashed lines are Hernquist profiles with scale radii different by a factor of two. Right: Scaled caustic mass profiles for the most 
massive quartile (ordered by M200) of CIRS clusters with the same model profiles. 



robust to the particular choice of interloper removal algo- 
rithms, while their virial masses have significant scatter 
relative to the CIRS values (primarily due to differences 
in estimates of r^oo)- The Lx — o v relation is expected to 
have the form Lx oc a p under simple theoretical models 
(Quintana & Melnick 1982). For clusters, most recent es- 
timates of the slope of this relation are consistent with a 
slope of 4. For galaxy groups, however, some authors find 
shallower slopes (Mahdavi et al. 2000; Xue & Wu 2000), 
while others find slopes comparable to clusters (Ponman 
et al. 1996; Mulchaey & Zabludoff 1998). In the Group 
Evolution Multiwavelength Survey (GEMS), Osmond & 
Ponman (2004) find a slope of 2.31±0.61 for groups, but 
the slope for a combined sample of groups and clusters 



is 4.55±0.25. In a similar vein, Mahdavi et al. (2000) 
find that the best-fit relation for a sample of groups and 
clusters is a broken power-law with a shallower slope for 
groups than for clusters. Unusual systems may affect 
the Lx — <J V relation: for instance, Helsdon et al. (2005) 
show that galaxy groups with apparently very small a p 
and large Lx may have larger <j p than earlier estimates, 
bringing these groups more in line with a steeper rela- 
tion. Khosroshahi et al. (2007) find that fossil groups 
follow a Lx — cr p relation with slope 2.74±0.45, shal- 
lower than clusters (they also note a tendency for fossil 
groups to have slightly larger Lx for a given a p than non- 
fossil groups). In the 400d-SDSS sample, we selected the 
groups based on their X-ray emission, whereas most pre- 
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Fig. 8. — Concentrations cioi=rioi/a of the best-fit NFW pro- 
files verus Mioi for the 400d-SDSS groups (solid hexagons) and the 
CIRS clusters (open circles). The dashed lines show the expected 
trend from (thick dashed line) simulations and a parametric model 
(thin dashed line Bullock et al. 200f). The solid and dash-dotted 
lines show the median value of cioi and log(cioi) for the CIRS 
clusters as a function of cluster mass. 




0.0001 0.001 0.01 0.1 1 10 

Lj, (10 44 h- z erg s"') 



Fig. 9. — Velocity dispersions at r2oo compared to X-ray lumi- 
nosities. The solid line is the bisector of the least squares fits for the 
CIRS clusters. The dashed lines show the a p — Lx relations from 
RASS-SDSS (Popesso et al. 2005). Symbol types are explained in 
Figure 1; large open circles highlight groups strongly affected by 
nearby structures. The two arrows indicate the locations of these 
data points for the larger clusters located close to two of these 
groups (A1541 near cll227+0858 and A1783 near cll343+5546). 

vious studies (Ponman et al. 1996; Mulchaey & Zablud- 
off 1998; Mahdavi et al. 2000; Osmond & Ponman 2004) 
begin with catalogs of optically selected groups and then 
analyze either pointed X-ray observations or RASS data. 

The 400d-SDSS X-ray luminosities are measured in the 
ROSAT band (0.1-2.4 keV) and corrected for Galactic 
absorption (see Burenin et al. 2007, for details). Fig- 
ure 9 shows the Lx — <J V relation for the CIRS clusters 
along with the best-fit Lx — cr p relation of the RASS- 
SDSS (Popesso et al. 2005) and CIRS clusters, the GEMS 



groups, and the GEMS groups+clusters fit (Osmond & 
Ponman 2004). The 400d-SDSS groups roughly follow 
the relations of the RASS-SDSS and CIRS samples, but 
the scatter is large. Squares, triangles, and crosses in Fig- 
ure 9 show systems from CIRS, Mahdavi et al. (2000), 
and Osmond & Ponman (2004) respectively. Note that 
the 400d observations are typically not deep enough to 
separate group emission from contributions from indi- 
vidual galaxy halos, so the X-ray luminosities are not 
necessarily directly comparable to those in some other 
studies (e.g., Mahdavi et al. 2000; Osmond & Ponman 
2004). However, because the 400d-SDSS groups lie close 
to RASSCALS and GEMS groups, this effect is probably 
small. 

The groups previously identified as undergoing merg- 
ers are highlighted with large circles in Figure 9. These 
groups are the most extreme outliers from the CIRS and 
RASS-SDSS scaling relations and the locus of points 
from earlier studies. This result can be easily under- 
stood as a result of confusing an Lx measured for a sub- 
cluster with the velocity dispersion a v reflecting the ve- 
locity dispersion of the parent cluster. Consistent with 
this explanation, the arrows in Figure 9 show where 
cll227+0858 (cll343+5546) would appear if we used the 
X-ray luminosity of the nearby, more luminous system 
A1541 (A1783); both points move closer to the CIRS 
Lx — o relation. Note that only one of the 400d-SDSS 
groups (ell 159+5531) is a likely fossil group (Vikhlinin 
et al. 1998), so the sample is not dominated by these 
unusual systems. Specifically, Voevodkin et al. (2009) 
search for fossil groups and clusters in the 400d survey; 
clll59+5531 is the only fossil group at z<0.1, and it only 
qualifies under a relaxed definition of fossil groups. 

Excluding the four groups contaminated by large-scale 
structure, the bisector of the least-squares fits has a 
slope of 3.4±1.6, intermediate between the shallow slopes 
found by some investigators for groups and the steeper 
slopes found for clusters (see Figure 9). Because the un- 
certainty in the slope for the 400d-SDSS groups is large, 
the slope is not significantly different from any of these 
relations. 

We test the consistency of the Lx — cr p relation for 
400d-SDSS groups with that of the CIRS clusters by fit- 
ting these relations for the CIRS clusters alone and for 
the CIRS clusters with the 400d-SDSS groups (exclud- 
ing those contaminated by large-scale structure) added 
to the sample. For the CIRS data alone, the slope of the 
bisector of least-squares fits is 3.7 ±0.6, and the slope for 
the combined CIRS and 400d-SDSS data is 4.1 ± 0.4. 

From the limited sample, we conclude that the X-ray- 
selected 400d-SDSS groups have an Lx — <J P relation con- 
sistent with previous studies of both larger-mass clus- 
ters (e.g., CIRS) and smaller-mass groups (e.g., GEMS). 
We find no evidence of large biases in this relation be- 
tween X-ray-selected and optically-selected systems. An 
important caveat to this conclusion is that some recent 
studies of clusters (Popesso et al. 2007) and groups (Ras- 
mussen et al. 2006) find that there may be some systems 
that arc significantly "X-ray undcrluminous," perhaps 
because these systems are dynamically young. 

Figure 10 shows Lx versus M 500 as estimated from the 
caustics for the 400d-SDSS groups. Figure 10 also shows 
the data from CIRS and the best-fit M 500 — Lx relations 
for CIRS (caustic masses) and RASS-SDSS (both virial 
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Fig. 10. — Left: Caustic masses at rsoo compared to X-ray luminosities. The solid line is the bisector of the ordinary least squares fits to 
the CIRS data. The red and blue lines show the M500 — Lx relations for RASS-SDSS (Popcsso ct al. 2005) for optical and X-ray masses 
respectively. Symbol types are explained in Figure 1; large open circles highlight groups strongly affected by nearby structures. The two 
arrows indicate the locations of these data points for the larger clusters located close to two of these groups (A1541 near cll227+0858 and 
A1783 near cl!343+5546). Right: Same for M200 - L x . 



masses and Tx-based masses Popesso et al. 2005). 

Figure 11 shows the M 2 oo — o> relation. The tight rela- 
tion indicates that the caustic masses are well correlated 
with velocity dispersion estimates. The good correlation 
is perhaps not surprising because both parameters de- 
pend on the galaxy velocity distribution. The best-fit 
slope for the CIRS clusters is M 20 o oc o-3.i8io.i9 with the 
uncertainty estimated from jackknife resampling. The 
dashed line in Figure 11 shows the virial mass-CT p relation 
for dark matter halos in simulations (Evrard et al. 2008, 
slope 3.035±0.023). The agreement between these rela- 
tions from CIRS and simulations further supports the use 
of caustics as a mass estimator (we compare the caustic 
masses to virial mass estimates in §3.5). The four merg- 
ing groups in 400d-SDSS represent the most extreme out- 
liers from the CIRS relation. The remaining 400d-SDSS 
groups approximately follow the CIRS relations, but we 
do not attempt to fit a relation for the 400d-SDSS groups. 

In the previous section we found that many of the 400d 
groups are either group-group mergers (cl0820+5645) 
or group-cluster mergers (cll033+5703, cll227+0858, 
cll343+5546). These groups represent the most extreme 
outliers in the Lx~<J p relation (Figure 9) and the Lx—M 
relations (Figure 10). This result reinforces the idea that 
a detailed understanding of the large-scale structure sur- 
rounding low-mass systems is critical to obtaining accu- 
rate scaling relations. Previous studies, (e.g., Helsdon & 
Ponman 2000; Mahdavi et al. 2000) usually exclude such 
systems for similar reasons. 

5. DISCUSSION 

We use the Fifth Data Release of the Sloan Digital Sky 
Survey to study galaxy groups and compare their prop- 
erties to more massive clusters. In particular,we study 
the presence of infall patterns around galaxy groups, use 
these patterns to measure mass profiles, and use these 
masses to compare the scaling relations of groups and 
clusters. This study extends our previous work on infall 
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Fig. 11. — Caustic masses at r2oo compared to velocity disper- 
sions within r2oo • The solid line is the bisector of the ordinary least 
squares fits for the CIRS clusters. Symbol types are explained in 
Figure 1; large open circles highlight groups strongly affected by 
nearby structures. The dashed line shows the virial mass-<x p rela- 
tion for dark matter halos in simulations (Evrard et al. 2008). 



patterns around X-ray clusters in SDSS (Rines & Diafe- 
rio 2006) to smaller X-ray fluxes more typical of groups. 
These systems typically have smaller X-ray luminosities 
and masses than the systems studied in previous inves- 
tigations. Despite their smaller masses (compared to 
CAIRNS and CIRS clusters), well-defined infall patterns 
arc present in most of these groups. Four of the 16 groups 
are contaminated by large-scale structure, preventing an 
accurate determination of their infall patterns. 

We use the infall patterns to compute mass profiles for 
the groups and compare them to model profiles. Group 
infall regions are well fit by NFW and Hcrnquist profiles 
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and poorly fit by singular isothermal spheres, similar to 
the results from CAIRNS and CIRS for clusters. Ob- 
served clusters and groups resemble those in simulations, 
and their mass profiles are well described by extrapola- 
tions of NFW or Hernquist models out to the turnaround 
radius. The shapes of the best-fit NFW group mass pro- 
files agree reasonably well with the predictions of sim- 
ulations; the average concentration is 0200^6.5, slightly 
smaller than the value for group size halos extrapolated 
from Bullock et al. (2001) and with similar scatter. These 
mass profiles test the shapes of dark matter haloes on a 
scale difficult to probe with weak lensing or any other 
mass estimator. 

There has been much discussion of whether galaxy 
groups follow different scaling relations than clusters, 
particularly with respect to the L x — & P relation (e.g., 
Mahdavi et al. 2000; Helsdon & Ponman 2000; Osmond 
& Ponman 2004) . Using the caustics to determine group 
membership, velocity dispersions, and mass profiles, the 
400d-SDSS groups generally follow the same scaling re- 
lations as clusters, but the small number of 400d-SDSS 
groups does not allow us to exclude results suggesting 
that groups have a shallower Lx — & p relation. 

Understanding scaling relations of clusters and groups 
is critical for future attempts to constrain dark energy 
with cluster abundance measurements. Detailed mass 
estimates are difficult to obtain for groups and low-mass 
clusters because their X-ray emission is typically faint 
and because they contain fewer galaxies than massive 
clusters. This work represents part of an effort to in- 
crease the sample of groups and low-mass clusters with 
detailed mass estimates. Future work is needed to mea- 
sure cluster scaling relations robustly across a wide mass 
range. In particular, deep high-resolution X-ray obser- 



vations of groups like those studied here could enable 
measurements of group masses and X-ray temperatures, 
while additional optical spectroscopy could improve mea- 
surements of velocity dispersions and dynamical masses. 
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